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The cosmology of general fourth order corrections to Einstein gravity is consid- 
ered, both for a homogeneous and isotropic background and for general tensor per- 
turbations. It is explicitly shown how the standard cosmological history can be 
(approximately) reproduced and under what condition the evolution of the tensor 
modes remain (approximately) unchanged. Requiring that the deviations from Gen- 
eral Relativity are small during inflation sharpens the current constraints on such 
corrections terms by some thirty orders of magnitude. Taking a more conservative 
approach and requiring only that cosmology be approximately that of GR during Big 
Bang Nucleosynthesis, the constraints are improved by 4 — 6 orders of magnitude. 

PACS numbers: 04.50.Kd, 98.80.Jk 



I. INTRODUCTION 

General Relativity (GR) is often considered to be an extremely well tested physical theory, 
however compared to the other fundamental constants (e.g. the electron charge, the speed 
of light etc.), the gravitational constants are rather poorly constrained. Whilst it is true 
that the predictions of GR have been verified to staggering accuracy (see for example (H), 
there is little experimental/ observational evidence to restrict corrections of GR that involve 
higher powers of curvature invariants. This is simply a consequence of the fact that GR has 
been directly tested in weak curvature regimes, with strong curvature tests (such as near 
black holes and neutrons stars) typically resulting in theoretical restrictions. The Einstein- 
Hilbert action for GR contains only terms that are, at most, second order in the derivatives 
of the metric (in the form of the Ricci scalar R) and a natural question to then ask is what 
restrictions are there on the presence of terms that contain higher order derivatives of the 
metric? 

Motivation for such a question can be found from two, broadly different, approaches. 
The first is phenomenological: are there corrections to the Einstein-Hilbert action that can 
better describe observed data? In particular there has been much effort in looking for devi- 
ations from GR that might explain the apparent presence of Dark Energy, Dark Matter and 
inflation (such as f{R) gravity j2l"12|, chameleon models [l3-15|, conformal gravity (Tg], 17 



MOND jl8l-l20[|. TeVeS [2lH26| etc.). The second approach is to consider theoretically moti- 
vated corrections to GR, often due to Quantum Gra vity (such as String/ Brane theory [27i - 



38], ADS-CFT |39|-|4l|. Loop Quantum Gravity l^lQ, Chern-Simons theory |52|, |53j etc.) 
or a Grand Uniflcation theory (such as Non- commutative Geometry j54-56||). Because there 
are only three curvature invariants that contain at most two derivatives of the metric, i?, the 
cosmological constant A and the Gauss-Bonnet combination (defined below), one is gener- 
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ally forced to consider correction terms containing 4*^^ and higher derivatives of the metric^, 
at least in the low energy, effective limit of the full theory. 

Theoretically there is, however, a serious difficultly with a theory that contains higher 
than second order derivatives of the metric; such a theory would contain 'ghosts'. These 
are modes that have a negative kinetic term and can result in super-luminal propagation 



and other instabilities |57H59j. Such difficulties become particularly acute when the theory 
is quantised, since the presence of ghosts leads to particles with negative energies or states 
with negative norm. This shows that any theory containing higher order corrections to GR 
cannot be a fundamental theory unless powerful non-perturbative effects come into play. 
Here we will consider a general 4*^ order correction to GR and take the point of view that 
this is an effective theory, approximating (some of) the corrections to GR that are produced 
by some complete underlying, non-perturbative theory. This precisely is the approach that 



is taken when renormalisation group techniques are applied to GR |6Ch62 . 

The only terms containing at most 4*^^ order derivatives of the metric that are constructed 
out of curvature invariants (other than the cosmological constant) are 

R'^'^^Rp.ufrr , R^^Rfiu , R^ , R , (1-1) 
where -R^^^p^ is the Riemann tensor, defined as 

pA* _|_ po pM pa pM {-\ n\ 

^ vpP ~ ^ pv,l3 ^Pv,p^^up^Pa ^ vl3^ pa 1 \^-'^) 

and we follow the notation of [63| Greek letters to denote space-time indices; is the usual 
Christoffel symbol; we denote partial derivatives as = R^ and we use the convention 
Rpv = R^^pu) with the signature ( — h ++). The Ricci scalar (last term in fll.ip ) is at most 
second order in derivatives of the metric, while the particular combination, R^ — ARp^i^R^" + 
R^'^p'^ Rp_yp^, called the Gauss-Bonnet combination, satisfies (in four dimensions) 



[R^ - 4:Rp^R^"' + R^'^P^Rp^p^) = total divergence . (1.3) 

Thus, assuming there are no boundary terms (which can be a significant complication, 
particularly when attempting to quantise the theory ^3]); "^^^ write any 4*^ derivative 
correction to the Lagrangian of GR as 

5 = - y dVv/^ [^i? - h(3R'' + haRp^R'^'^'j + Sm.tter , (1.4) 

where for simplicity we have neglected the cosmological constant. Note that hig her order 



terms (e.g. R^) can lead to 4*"^ derivative terms in the equation of motion 6J|, however 
here we will restrict our attention to terms that are of similar order in the Lagrangian, 
assuming that Eq. (11. 4p is a low curvature expansion of some underlying theory. We define 
= 327rG where G is Newton's constant and we use units in which [c] = 1. a, P and 7 are 
dimensionless couplings and the presence of h is due only to the dimensions (i.e. the theory 
is entirely classical). In general the coefficients a and /3 are arbitrary constants, however if 
these corrections are due to some underlying Quantum Gravity theory, then they would be 
expected to be of order one, so that the quantum gravity scale is of the order of the Planck 
scale. The value of a and /3 set the scale at which significant corrections to the GR occur 
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and in the following we will restrict these values via cosmological considerations. In all of 
the following we set [^] = 1. 

By varying this action with respect to the metric g^^, as usual, one finds the field equa- 
tions (ssj . 

H^, ^ {a- 2/3) R^.^, - ai?'^,^^^ - (| - 2/3) g^^^R^ + 2aR^^R^^^^ 

-2(3RR^ - [aR'^Rpx - m') + l^^'G^u = -\t^u , (1.5) 

where R.^^ = V^i? is the covariant derivative of R and T^j^ is the energy momentum tensor 
given by varying iSmatter with respect to g^^. Perturbations of this theory around a fiat 
background, show that the theory contains, in addition to the massless graviton, a massive 
spin-2 field (of negative energy) which is a ghost field and a (positive energy) massive scalar 
field |58|. 

One particular example of a theory that predicts 4**^ order corrections to GR is Non- 
commutative Geometry jsH. Here the asymptotic expansion of this (non-commutative) 
geometric theory produces the entire standard model coupled to a gravitational action of 
the form (56| . 

5ncg = - y d^xv^ (l^^ + ^ncgC'^'^C^,^^ , (1.6) 

where C^up-y is the Weyl tensor^ . The consequences of this modification for cosmology (65i 
and astrophysics j63, 68 1 have been considered and in particular in [66^ it was shown that 
background cosmology is unaffected by the presence of such a correction. Using (11. 3p to 
write action this in the form of (11.41) gives. 



<Sncg = - d^xv^ 



1 A 1 2 

R + 2q;ncg I RpuR'^" — -R 



16nG 



[1.7) 



In particular then corrections of this form satisfy a = 3/3 and in the following we will show 
that only this combination leaves the background dynamics unaffected. 

The question we will address is what restrictions on a and /3 can be deduced from the 
cosmology given by (II. 4p and how they compare to the current constraints. Existing restric- 
tions on a and /3 are in fact very mild, requiring only that these coefficients be less than 
0(10''^). The restrictions come from solar system tests, notable the perihelion precession 
of Mercury j58| and gravitational wave production in binary systems [67|]. There are also 
laboratory scale tests of the inverse square law, which restrict parameters of this form signif- 
icantly more, a < O (10^°) [69.]. However here we will focus on the cosmological implications 
of such corrections and so compare to the other large scale constraints. 

The reason that these restrictions are so weak is that corrections to General Relativity 
occur at higher orders in curvature and hence are highly suppressed in weak curvature 
systems. In order for the effects of such corrections to become more significant, one needs 
to consider a system that contains either strong curvature (such as black holes or the early 
universe) or long evolution times (so as to allow the small deviations from General Relativity 
to accumulate). Here we consider the latter case, by examining both the dynamics of the 



^ It should be mentioned that this Non-commutative Geometry theory is formulated in a Euclidean signature 
and is entirely classical. 
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background cosmology and the evolution of tensor mode perturbations over cosmological 
time scales. 

Ideally one would want to consider the evolution of scalar mode perturbations, since these 
can be directly related to the wealth of observational data coming from large scale galaxy 
surveys, weak lensing maps and the Cosmic Microwave Background (CMB). By contrast, 
cosmological tensor mode perturbations, although generically predicted by inflationary mod- 
els, have yet to be observed. Despite this, tensor mode perturbations have the advantage 
of being technically much simpler to calculate than scalar modes (evolving according to a 
single equation, rather than four coupled equations). In addition, since tensor perturbations 
do not couple to the matter content of the universe other than through the background 
evolution (at least for matter with vanishing anisotropic stress) they are a direct probe of 
the underlying gravity theory and are unaffected by any possible modifications to the matter 
action. 

In Section ini we derive the dynamics of the background (i.e. homogeneous and isotropic) 
cosmology, discussing in particular the solutions during inflation and during the radiation 
and matter dominated eras. Section IIIII derives the general evolution equation for tensor 
mode perturbations for the action given in (11. 4p . while in Section HVl we consider the evo- 
lution of such perturbations against specific background cosmologies. Using these evolution 
equations we derive constraints on a and /3 in Section |V] and summarise and conclude in 
Section IVIl 



II. BACKGROUND COSMOLOGY 

In order to examine the background cosmology of this theory, we need to calculate the 
modified Friedmann and Raychaudhuri equations, which are given respectively by, = 
^T°Q and H\ = -^T'^i- We want to consider a Friedmann- Robertson- Walker universe, 
whose metric is given in block diagonal form by, 

where i, j ■■■ = 1,2,3 will be used to denote spatial components; is the three dimensional 
zero vector; a{r]) is the scale factor, with rj conformal time, and the three dimensional spatial 
metric is given by, 

7.. = -7-^ 72 > (2-2) 

where 6ij denotes the three dimensional Kronecker delta and K the (scaled) spatial curvature 
{K = 1 has spherical spatial topology and corresponds to a closed universe, K = Ois spatially 
flat and K = —1 is spatially hyperbolic). We shall denote background quantities with a bar. 
For example, the full metric is given by, 

g^^ {r], x) = g^^ {r]) + 6g^„ {r], x) , (2.3) 

where, as expected from homogeneity and isotropy, the background quantity (in this case the 
metric) depends only on conformal time, whilst the perturbed quantity (here Sg^i,) depends 
on both spatial and temporal coordinates. 

With this metric one can immediately calculate the components of the curvature tensors. 
For later convenience we list some of the non-zero components the of Riemann and Ricci 
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tensors and the Ricci scalar, with the remaining non-zero components being given directly 
from the symmetries of the tensors. 

iOj ~ lij ' OOj ~ ! 

- 3 

R'j = {n' + 2 (-H^ + K)) 6] , R = ^ i'H' + n^ + K) , (2.4) 

where a dash means differentiation with respect to conformal time rj and "H = a' /a is the 
con/orma/ Hubble parameter (see for example f70|). 

As outlined in Appendix A, one can use this metric to find, 



H 



2 {a- 3/3) 





0—2 -2 



~2H in" + 2HH') + Sn^ + 271^ K + {H'f - 



+ (^' + ^) = ^«V, (2.5) 

+ {W + 'H^ + K)= TTT^a' (p - 3P) , (2.6) 

where we used the fact that by homogeneity and isotropy the most general background 
energy momentum tensor is that of a perfect fluid, 

(2,7) 

with p and P the energy density and pressure of the fluid. 

Notice that, as expected, the a = P = case i.e. GR, gives the usual Friedmann and 
Raychaudhuri equations, as does a = 3(3. Recall that latter case corresponds to the NCG 
motivated correction to GR (see Section |T]), which vanishes for a FRW universe [66!]. 

Unfortunately there does not appear to be a general, analytic solution to (12. 5 p for a 
fluid satisfying P = up, with u a constant, however progress can be made by checking 
the consistency of certain, important, cosmo logical solutions. In particular we will consider 
radiation and matter dominated universes, whose matter components have a; = 1/3 and 
a; = respectively and also a universe dominated by a slowly rolling scalar field'^, for which 
u = —1. In [tiI the conformal properties of (11. 4p were used to derive general consequences of 
these corrections, however here we explicitly calculate the background evolution, to facilitate 
the inclusion of tensor perturbations in Section IIIII 

One would like to be able to specify the form of the matter content (i.e. u) and deduce the 
corresponding cosmological history a{ri), however this leads to rather involved expressions 
and a more transparent approach is to specify the form of the cosmological expansion a (77) 



^ Recall that the potential energy and pressure associated with a scalar field evolving in a potential V (0) 
are — 0^/2 + V {4>) and = — V {(p) respectively. Thus if the field is slowly rolling i.e. if 

V {(p) >• 0^/2 then the pressure is given by ~ —p^ 
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that matches the GR expectation and evaluate the properties of the matter components that 
support this expansion. The properties of these matter fluids will, in general, deviate from 
their standard values and this deviation measures the significance of the correction terms 
on the cosmological evolution. It is important to note here that the non-linear nature of 
Eq. (12. 5p and Eq. (12. 6 p means that finding such deviations to be small, does not immediately 
imply that the solutions, a{ri), with the correct matter content would be close to those of 
standard GR. The approach used here is however a necessary (but not sufficient) condition 
for the cosmological evolution to approximate that of GR and should be supplemented by 
direct numerical evaluations of Eq. (12. 5p and Eq. (12. 6 p . 



A. De Sitter expansion 

The scale factor for an exponentially expanding universe (in conformal time) is given by 

«W) = ^— > (2-8) 

where Oinf is the normalisation of the scale factor at the beginning of inflation, taken to 
occur at (conformal) time r/inf. Using this ansatz in (12. 5p and (12. 6 p we find. 



12 (a - 3/3) 

(ainf'7inf)^ 



p{a) 



(Oinf^inf)^ V V a / 

p („)(!_ 3.) = ^Ih^mK^^^J^^Kr-^)'] . (2.9) 

(ainf??inf) a? (ainf?7inf) V \ a J J 

In a spatially fiat universe {K = 0), we then have 

p{a) = -^ ^1-87 72] , (2-10) 



ainf'7inf) V (ainf?7inf)' 

and 

^ = ^1-4(1-8 ^ ) , (2.11) 

-3 \ V (ainf?7inf) / / 

where 

A = ^. (2,12) 

Recall that both k and cosmic time t = at] have units of length (or equivalently M^^), hence, 
on dimensional grounds, one expects the key dimensionless parameter to be A {ari)~'^. 

From (I2.10p we see that exponential expansion is given by a fluid with a constant energy 
density just as for GR, however for a given energy density the Hubble rate, Ti, will depend on 
the parameters of the theory [a, /3 and 7). In order to maintain the exponential expansion 
the equation of state parameter of the fluid, u differs from the GR case (a; = — 1). Although 
it is a constant, it is important to note that this is only true for exactly exponential expansion 
with = 0, in general one would require a varying u. The equation of state that is required 



-1.5 -1 -0.5 0.5 1 1.5 

^(ainf?7inf)~^ 

FIG. 1: The equation of state that is required to maintain an exponential expansion depends on the 
coefficients of the corrections to GR, A (ai^^'qi^^)^'^ . Also plotted for comparison are the equations 
of state for a cosmological constant cj = —1, pressureless matter u = and radiation uj = 1/3. 



depends on the a, (5 and 7 and is plotted in Fig. ([I]), notice in particular for SA = (oinfr^inf) 
the equation of state that produces exponential expansion diverges, essentially meaning that 
it is not possible to have such expansion for these for a fluid with a finite pressure. However 
for 8^4 = (ainf77inf)^, f l2.10p gives p (a) = 0, thus it is not surprising that the pressure becomes 
ill defined. For this specific choice of parameters, the vacuum solution of the theory is an 
exponentially expanding cosmology. This divergence will appear again (for the a = case) 
in Section llllt when perturbations are considered. 

From Fig. we see that for < (ainf77inf)^^ < 1 a phantom fluid is required to 
produce the expansion i.e. u < —1, while for (ainfr^mf)^^ > 1 a fluid with u > 1/3 is 
necessary. Finally for A < (i.e. a — 3/3 < 0) exponential expansion can be achieved by 
fluids with — 1 < w < 1/3. In particular, only A = {or equivalently a — 3/3 = 0) will 
produce an exponential expansion in the presence of a slowly rolling scalar field (for which 
u ~ —1). The A < and 8A (ainf?7inf)~^ > 1 cases could have the desired expansion driven 
by a scalar field if the kinetic and potential energies of the field were correctly matched, to 
give (12.111) . However even if this does not result in excessive fine-tuning (over that required 
by standard infiation), it would likely be rather difficult to ensure that u has the correct 
value over a long period of expansion. 

Thus a period of exponential expansion, within general fourth derivative gravity theories 
is unlikely to be sourced by a simple scalar field, unless A -C 1 (this is the case for example 
in j72| and j73|). However the presence of non- minimal couplings between the scalar field 
and the curvature terms can have a significant effect even in the a — 3/3 = case 



so one cannot rule out the possibility. Indeed, the motivation for studying such corrections 
to GR requires that we consider a — 3/3 be small (and hence A to be small), in order to avoid 
super-luminal propagation (and other ghost effects, see Section [T|). With this in mind it is 
possible that standard slow roll inflation would result in almost exponential expansion i.e. 
that the deviations from the expected GR behaviour are small. Determining the observa- 
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tional consistency of such an approximation would require both a general solution to (I2.5P 
in the presence of a scalar field and, eventually, details of the evolution of perturbations 
against this background. 

If we want to have inflation driven by a scalar field in the standard way, the correction 
terms need to satisfy, \A\ (amf^^mf)^^ ^ 1- Since this is a constant, it would imply that also 
at the end of the inflationary epoch |y4| (ainf^^inf)"^ ^ 1 and hence the radiation era would 
begin with this dimensionless ratio being small. 

B. Radiation era 

In standard GR an era dominated by a single fluid with u = 1/3 (i.e. radiation) produces 
a [r]) = aT-rj/rir-, where is the scale factor at r/r, which is the conformal time at which the 
radiation era begins. Using this form of the scale factor in (12.51) . we find that 

/ X — V- 3/3) ( ^j^fa^r 
^^"^ = [^^^""[-^ 

(2.13) 





and 



48 (a- 3/3) M , 
p (a) (1 - 3a;) = ^ ^ ^ + ^—^ . (2.14) 

Thus, if we restrict our attention to the K = ^ case, we find that an expansion of a [j]) = 
arTj/rjr is produced by a fluid with a; = 1/3 and 

where A is given in f l2.12p . Such a fluid has the equation of state of radiation, however the 
required energy density would be larger (for a — 3/3 > 0) than that of radiation and would 
scale differently. Such corrections to the energy density would, however, rapidly become 
small as the universe expands and 

8 < 1 • 2.16 

As discussed in Section f ill Ap . if the radiation era is preceded by an epoch of inflation that 
is driven by a scalar field, then 8|A| <C (ctr^r)^ at the beginning of the radiation era. This 
would ensure that the expected expansion would approximately be achieved by standard 
radiation even early in the radiation era. 



C. Matter era 



Repeating the above procedure for matter era type expansion, i.e., 0(77) = a^rf/rj^ 
where as in the previous two cases, we have defined the time the matter era began as rj^ 



2 

m' 
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and the scale factor at this time to be a^^. One then finds (again for the = case) that 
the required energy density is, 

Thus, just as in the radiation case, the fiuid would have to have a larger (for a — 3/3 > 0) 
energy density than standard pressureless matter. These corrections would rapidly become 
very small as the universe expands, however, unlike the radiation era case where the usual 
equation of state was found, here we require, 

- = ^(l + 344?^V' . (2.18) 



Therefore an exactly pressureless fiuid will not produce an expansion consistent with pre- 
dictions of GR. However the pressure that is required to match the GR expectations is 
vanishingly small for 

^ « 1 . (2.19) 

Thus for large scales (or late times) the background cosmological evolution will be ap- 
proximately the same as that expected from GR, provided 

\A\ (ainfr/inf)"' < 1 . (2.20) 

Considering for the moment only the radiation and matter eras, all that is required is 
that \A\ (ar?7r)~^ be sufficiently small, so as to allow fl2.15p to be approximately correct for 
standard radiation. As the universe expands the evolution will then quickly approach the 
standard cosmology of GR. For \A\ (ainf'7inf)~^ ^ 1 this is satisfied and hence the standard 
background dynamics can be (approximately) reproduced from a slowly rolling scalar field, 
with the corrections to the GR expectations remaining small throughout the history of the 
universe. This agrees with the conclusions of (7l| . 

We have considered only the fiat {K = 0) case, the presence of a non-zero curvature would 
introduce additional corrections to the standard GR evolution equations. However it is clear 
from (12. 9p . fl2.13p and f l2.14p (and similarly for the analogous equations for the matter era) 
that the consequences of i^' 7^ rapidly become small in a large scale, expanding universe. 
This is to be expected since the energy density associated with the spatial curvature of an 
expanding universe decreases with time, hence its consequences will diminish. 



III. TENSOR MODE PERTURBATIONS 

In the previous sections we derived the evolution of the background cosmology for a 
general 4**^ order correction to GR. Given the prevalence of ghost modes in such theories, 
it is important to examine the evolution of perturbations on these background cosmologies. 
Here we derive the equations governing the dynamics of cosmological gravitational waves. 



See [75[ for an alternative derivation in a more general setting and [76[ for an derivation 
similar to the one presented here. Also [t^] for the generation and evolution of tensor modes 
during infiation for a specific f{R) model, which includes a detailed analysis of the properties 
of the resulting spectrum. 
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A general perturbation to the metric, fl2.ip . can be decomposed into scalar, vector and 
tensor components (see for example [t^) and here we will be concerned only with the tensor 
modes, which are transverse and traceless and given by the perturbed metric, 

'^fl'oo = ^g^i = , 5gij = 2a^ (?7) Eij , (3.1) 

where the perturbation Eij is transverse {D^Eij = 0) and traceless {g'^^Eij = 0), the non- 
transverse and non-traceless parts giving contributions to the vector and scalar perturbations 
respectively. Using this perturbed metric one finds the perturbations to the components of 
the Riemann tensor, 

^R^ioj = ~ i^ij + 'HE-j + 271' Eij) , SR^ij^. = E[j\f^ — E-^i^- , 
^^^00, = {E")\ + -H {E'Y^ , 5R\^, = - (E')\|, + {E'Y^^, 
SR',ok = - iE')J' - (E')\|, , 

-(e'i -E'-, +E\- -E' r +E. -E.^\ 

\^ j\nm ^ j\mn ' ^ n\jm ^ m\jn ' jm \n jn \m 

- -H {61 (E'),. - {E')^^ - {E')\ 7,™ + {E')\ 7,n) , (3.2) 

where Eij\k = D^Eij is the spatial covariant derivative of Eij. Similarly one can calculate 
the perturbations to the remaining tensor components given in 02.41) . the only non-zero 
expression being, 

SR^^ = 1 (- {E"y^ - 2V, {E')\ + {A-K) E^) , (3.3) 

where A = D^Di is the spatial Laplacian. Note also that, because Eij is traceless, we have 
5g\ = 0, and hence 6G^j = SR^j. 

With this one can calculate the perturbations to fll.5p to find, 



= (a - 2/3) 5 {R^-.^J) - aS (i?'^,^^^,) - (| - 2/?) g^J {R%) + 2a {SRY' R 



pX jyfi 

puX 



+2aR'^' {6Rr. + {^K-' - 2m) mY. = -i^sn ■ (3.4) 



Of the various terms appearing in (13. 4p those not involving covariant derivatives are the 
easiest to calculate. Using (13. 2 p and (13. 3p . one immediately finds. 



(SR)'^ R\j, = J2^^±^ hE'y^ + 2n {E'y^ - {2n' + + a + 2k) , (3.5) 

and 



R'\^RY„x = ^,{-sn'iE"y^-2n'n{E'y. 

+6 {H^ + K) H {E'y^ -[W + 2 [H^ + K)) [2H^ + A) E'^ , (3.6) 
while since only tensor perturbations are being considered the remaining components vanish. 
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Calculating the terms in (13. 4p that involve covariant derivatives is rather more involved 
however it is easily achieved by expanding out the covariant derivatives in terms of the 
metric and its derivatives. Here we quote the results, with the (tedious) calculation given 
in Appendix B. 



(jj 

+ (~ ^^"^^^ ~ ^^'^^^ + i^-K) i?',)) , (3.7) 

where in the last equation, the differential operator V is defined as 

V^(^-^^-2nd, + 2'H' + A^ . (3.8) 

Using the expressions given in (13. 7p one readily finds that the evolution equation for 
tensor mode perturbations, (13.41) . becomes 



-a {E""y^ + (ai + 2aA) {E")\ + {E')\ 



+ [(A - K) (ag - aA) + - 2a {W + 2%^ + 2K) A] E' ^ 
where we have defined the rj dependent functions. 



(3.9) 



ai = -AW {2a + 3/3) + {a - 3/3) - 3K {a + 4/3) - 



7a 



2 



K 



2 



as = 4H" {a - 3/3) + SaU^ + QaUK - mW {a + 6/3) - 



270^^ 



as = -2'H'{a-6f3) + 12f3W + 2K{a + 6f3) + 



7a 



2 



K 



2 



a^ = Qan'^K + AaKn' + QaK^ . (3.10) 

Notice in particular that the degeneracy present in the background equations between 
a and /3 is now broken and the particular combination a — SP play no special role in the 
perturbation equations. This agrees with the results of [66, which showed that deviations 



from standard cosmology for NCG inspired corrections to GR (for which a = 3/3), are 
important for perturbations in general and gravitational waves in particular. Indeed, by 
restricting our attention to fiat Minkowski space-time, i.e., "H = and K = 0, and taking 
a = 3/3, (13. 9p reduces to the equation for gravitational radiation derived in j67j . 

IV. TENSOR PERTURBATIONS AGAINST DIFFERENT COSMOLOGICAL 

BACKGROUNDS 

The general evolution equation given in (13. 9 p can be simplified by considering specific 
forms of background evolution. For example, consider a {rj) = a\ri^"r]'^ for some constants z/. 
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a-i and 77;. Such power law behaviour is typical in cosmologies dominated by a single matter 
fluid in GR e.g. v = 2 corresponds to the matter dominated era, while v = 1 gives the 
radiation dominated era, with a\ and rji being the scale factor and conformal time at the 
beginning of these eras respectively (see Section HT])- One can easily extend this to include 
several matter fields and hence allow for (for example) the transition from the radiation to 
matter dominated eras, however for simplicity we focus on single fluid systems. Here we 
have kept a and /3 general, however it must be remembered that if (for example) a [rj) oc r) 
is to represent the expansion during a radiation dominated era, then we must satisfy fl2.16p . 

Decomposing all spatial functions into eigenf unctions, Qk of the spatial Laplacian i.e. 
AQk = — IkpQk ^ so that in particular AE^ = — |kpi?^, one finds. 



6W 



powerlaw 



|k| 

a; 



X 



—a- 



dx^ 



X 



+ I + ^ - /C 

\Xi 



+ I 6ai^^ + ^-2uK: 



X"^ \Xi 



X 



Xi 



K 




(4.1) 



where 



/3i = 2a + 3 (a + 4/3) 



K 



/32 
/36 



4z/(2a + 3/3 + z/(a-3/3)) , 
Sz/^a + 4z/2 (a + 6/3) + 8z/ (a - 3/3) 



4 (a - 3/3) 



K 



+ 4a 



6z/ (z/ (a - 2/3) - a + 2/3) 
K 



K 

w 

K 



a 



2v [a - 6/3) - 12Z/2/3 
+ 2ai/(2z/- 1) , 



-3 (a - 4/3) 



(4.2) 



and we have changed variables to x = ?7|k| and defined the (dimensionless) variables x^ = 
|k|?7i and 



(4.3) 



Note that if /C S> a,/3, then (14. ip approximates the usual GR evolution equation for tensor 
modes (e.g. [to!]), as it should. In particular then the evolution of extremely large scale modes 
(i.e. modes for which |k| <^ 1) is precisely that of GR. Of course, here one has to be careful, 
since the GR limit of Eq. (14.1 p reduces the order of the perturbation equation and hence can 



^ For the spatially flat, K — case this is just the usual Fourier decomposition. For K ^ 1 the eigenvalues 
take discrete values, |kp = I {I + 2) for I E Z and for K = —1 the eigenvalues are bounded below by 
|kp > 1 [70]. 
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introduce strong instabilities. One can, however check that 2/ the solution matches the GR 
solution initially, then it will continue to do so throughout the evolution. Furthermore, one 
can estimate the magnitude of the effect of small deviations to the GR solution, by evaluating 
the higher derivative terms appearing in Eq. (14. ip for the (approximate) GR solutions. Doing 
this one finds that the corrections seem to remain small during the evolution of the modes, 
however such evidence does not remove the possibility of their being a growing mode to the 
corrections to the GR solutions, even in the K, ^ a, P limit. However, if such an instability 
were present, it would likely involve a significant growth in the amplitude of and hence 
would result in even stronger constraints than those presented below. 

Of particular interest is the evolution of tensor modes in a flat [K = 0) universe, during 
the radiation era. Recall for such a radiation dominated universe we have u ^ 1 (see the 
discussion below f l2.15p ). One immediately sees that for v = 1 the coefficients given in (14. 2 p 
are independent of /3. Hence the propagation of tensor modes during the radiation era are 
unaffected by the presence of B? corrections, for a fiat universe, regardless of the strength of 
the coupling /3. Assuming that the universe is exactly flat, the evolution during the radiation 
era {v = 1) for a general a is given by. 



5H\ = a{'-^\ - 




(4.4) 



Hence, the evolution will be well approximated by that of GR when 



X 



a \x 



> 1 



(4.5) 



which is true at late times, i.e., large x, (recall however the comments below Eq. (14. 3p ). 

If we consider the evolution of perturbations to a fiat inflationary era i.e. with K = Q 
and V = —1, and consider only corrections i.e. a = 0, we find that (14. ip factorises into 



6W 



-1, 0=0, K=0 



24/3 



d^E% 2 dE\ 

X ax ■' 



da;^ 



(4.6) 



The differential equation in parenthesis is precisely that of General Relativity, thus, in 
the absence of a source, the evolution of tensor mode perturbations during inflation (with 
i^' = 0), is unaffected by the presence of corrections (provided there are no R^^R^^, 
corrections). Note however that although the evolution of these modes is independent of 
such corrections, their generation will not be. As can be seen from (14. 6 p the effect of the R^ 
correction is to introduce an effective, Newtonian constant 



7 



— )■ tt; 



7 



+ 24 



(4.7) 



In particular there is a critical value at 24/3/?^ = —7 (ai?7i)^ at which the evolution becomes 
ill defined. This is precisely the same divergence found for the background evolution in Sec- 
tion pi Ap . The vanishing of the prefactor in (14. 6 p makes the evolution arbitrarily sensitive 
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to any deviations for this equation of motion (i.e. the evolution is unstable), which is to be 
expected for perturbations about such a pathological background cosmology. 

We mention in passing that for the specific ratio a = 24/3 and K = 0, (14. ip can be ex- 
actly solved. The resulting expression, involving hypergeometric functions, is not, however, 
particularly illuminating. 



V. CONSTRAINING MODIFICATIONS TO GR 

As mentioned in Section [H local tests of GR, particularly from the perihelion precession of 
Mercury [s^ and the production of gravitational waves from binaries j68| , provide only very 
weak constraints, typically of the order of (a — 3/3) < 10™. Cosmology has the advantage 
over such test, both because it is possible to reach higher curvature regimes in the early 
universe and also because we naturally have very long evolution times, during which even 
small deviations from GR can become significant. As we found in Section [Tll in order for 
a slowly rolling scalar field to have sustained the exponential expansion of inflation, in the 
usual way, we require 

\a - 3I3\\ / 1 ^ ^ 



, . , «1 ' (5-1) 

where recall ainf and ?7inf are the scale factor and conformal time at the beginning of inflation. 
Note that r] is conformal time, so ainf?7inf = ^inf, where tinf is the cosmic time at the beginning 
of inflation. Taking inflation to have occurred t ^ to x 10~^° ~ 10~^^sec [78], where to is the 
current cosmic time (taken to be approximately lO^^sec), we see that (15. ip becomes. 



\a-3(3\ 



or, taking 7K ^ = (IGvrG) ^ and using the experimental value of Newton's constant {G ~ 
2.88 X 10~®''sec^ in these units), 

|a - 3/3| < 6.7 X 10^° . (5.3) 

While this constraint may still seem rather weak it is an improvement by some thirty or- 
ders of magnitude over the previous best astrophysical constraints (and is twenty orders of 
magnitude better than laboratory tests). 

One may be concerned that this constraint is produced by assuming that inflation is 
produced by a slowly rolling scalar field. Maybe there is some other mechanism for producing 
the exponential expansion required to solve the horizon, curvature and monopole problems 
of the standard big-ban g th eory? Indeed, there are examples of higher derivative theories 



driving inflation directly [65|, |79| . Or perhaps there is an alternative to inflation that does not 
require such exponential expansion? Certainly, since we have yet to derive the consequences 
of this general 4*^^ order theory on the growth (or otherwise) of scalar perturbations, it 
is premature to suggest that such an expansion phase would produce the scale invariant 
perturbation spectrum required by the CMB. 

To avoid these difficulties let us take the more conservative point of view and assume 
that the dynamics of the universe closely match those expected from GR only after some 
time, tc, in the radiation era. Certainly the success of the CMB tells us that GR is a good 
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approximation significantly before decoupling (i.e. before the last scattering surface). From 
f l2.16p we then have 

(5.4) 

where tr is the cosmic time at the beginning of the radiation era and ar and Oc are the scale 
factors at cosmic times and tc respectively. Being extremely conservative, let us assume 
that the universe was radiation dominated at energy scales below electro-weak unification 
(^Ew ~ 10~^^sec) tj- ~ 10~^^to^- Further let us assume that the universe is well approximated 
by that of GR at the scales on which Big Bang Nucleosynthesis takes place, tc ~ 10~^^to [ISj. 
This assumptions gives 

<10-^°t^ (5.5) 



7K 

or, putting in the current age of the universe and Newton's constant, |q; — 3/3| ^ 10™. Thus 
even with these extremely conservative estimates, the constraint is the same as that from 
gravitational waves and solar system tests. 

While the estimates above can provide a significant improvement on existing constraints, 
in the background equations there is a degeneracy in the parameters, in that we can restrict 
only the combination a — 3/3. In order to place independent constraints on each parameter, 
we need to consider perturbations around the background cosmology. Using the results of 
Section HVl we see that, during the radiation era, the evolution of the cosmological gravita- 
tional waves will deviate significantly from the predictions of GR unless. 



> 1 . (5.6) 



As a specific example, consider a mode that crosses the Horizon at matter-radiation equal- 
ity^, which occurs at teq ~ lO^^sec, and again conservatively consider the beginning of the 
radiation dominated era to be at the electro-weak scale. One then finds that (15.61) reduces 
to, 

a < 6.7 X 10^° , (5.7) 

again comparable to the current best constraints on such parameters. As is clear from 
(15. 6p . this constraint can be improved in two ways; by decreasing x/xi i.e. considering the 
modes close to the beginning of the radiation era or by decreasing /C. Recall from (14. 3 p 
that /C oc |k|~^, thus the constraint can be improved by considering smaller wavelength 
gravitational waves (i.e. larger |k|). The fact that these deviations from GR become more 
significant at small scales might have been expected because of the fact that these modes 
are ghosts • 

Taking the point of view that this theory is only an effective theory, good down to some 
scale, after which some more complete theory will dictate the evolution of these modes, we 



^ Here we are estimating the cosmic time at which this occurs by the GR value. Since we are assuming 
that GR is not a good approximation until we reach times beyond > t,-, this is not strictly consistent. 
However we are only estimating orders of magnitude, for which this should be sufficient. 

^ Close to matter-radiation equality, one should use a two fluid approximation, however we are concerned 
only with rough estimates and hence neglect this complication. 
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require only that the deviation for GR be not significant (i.e. that the ghosts are not appar- 
ent) on observational scales. The use of the word observational here deserves some comment, 
since we have yet to observe cosmological gravitational waves at any scale! However the de- 
viations from GR will typically make the tensor perturbations grow (as they represent an 
instability). If such growth had occurred and been significant, we would have observed their 
contribution to the galaxy power spectrum and the CMB. A quantitative constraint would 
require an in-depth analysis of the exact evolution given in f l4.4p and an understanding of the 
rate of the growth of the modes that can be expected. However in order to get an estimate 
on the order of magnitude of the constraint, we will simply require that the deviations from 
GR be small at the smallest observable scales at which such cosmological tensor modes could 
have been observed i.e. for the first modes that entered the horizon that have been observed 
in the galaxy power spectra. 



Data from Lyman alpha forests [8l|, |82| provide us with a detailed power spectrum 



of (scalar) density perturbations that matches the expectations of GR down to scales 
|k| ~ (10^ — 10^) |k|equaiity5 whcrc |k|equaiity IS the scale at matter radiation equality (which 
produces the characteristic peak in the power spectrum). Thus at such scales the constraint 
is improved by a factor of 10^ — 10^ over (15. 7p . This improvement would be further increased 
by considering smaller scales or, of course, by direct measurement of the cosmological grav- 
itational wave background. 



VI. CONCLUSIONS 

We have shown that it is possible to approximately reproduce the key epochs of standard 
cosmology within a generally modified 4**^ order gravity theory. More precisely we demon- 
strated that only a correction to Einstein-Hilbert gravity of the form C^'^^'^^C^ufri introduces 
no deviation from standard background cosmology, while any other correction can approxi- 
mate the standard background evolution, provided it is initially 'close' to such a correction 
(see Section ini for a more precise statement). Furthermore by requiring that the deviations 
from GR be small at certain key eras (such as during inflation and Big Bang Nucleosyn- 
thesis) we can quantify how 'close' general corrections need to be in order to be consistent. 
Even taking the conservative point of view that radiation domination began only after the 
electroweak scale, such constraints are comparable to the best current restrictions on the 
theory. Further requiring that inflation be sourced by a fluid with a; ~ — 1 allowed us to im- 
prove this constraint by thirty orders of magnitude] The appearance of such large constraints 
is not uncommon within cosmology. Consider for example the spatial curvature within non- 
inflationary cosmology. If the energy density associated with such curvature vanishes to one 
part in 10 today, it would have to vanish to one part in 10^° close to the Grand Uniflcation 
scale. Indeed this was one of the original motivations for inflation. Similarly requiring that 
the higher order curvature corrections be small today is a much weaker constraint than 
requiring they be small in the past. The earlier we can examine the consequences of such 
corrections (such as during the radiation era or inflation) the more the constraints can be 
improved. 

At the background level we are able only to constrain the parameter a — 3/3 (a = 3/3 corre- 
sponds to the C^^f^C^up^ correction of Non-commutative Geometry) and not the individual 
coefficients. In order to break this degeneracy we consider the evolution of perturbations 
against the homogeneous and isotropic background. For comparison with observations, one 
would like to examine the scalar perturbations of such a theory, which can be directly related 
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to the galaxy power spectrum, weak lensing surveys, and the CMB. However the equations 
governing the evolution of scalar perturbation are extremely involved in a general 4*^^ order 
gravity theory, making their use difficult. Instead we focused on the evolution of tensor 
modes (cosmological gravitational waves), which are technically simpler to calculate. While 
they have yet to be observed, they have the additional advantage of not coupling directly to 
matter (in the absence of anisotropic stress). Thus they are independent of any corrections 
to the matter action and are direct probes of the underlying geometry of the problem. 

Without direct observational evidence of such tensor mode perturbations, we considered 
the types of gravitational theories for which the evolution is approximately that of GR. 
The fact that such theories will typically lead to instabilities of perturbations, suggests that 
deviations from GR will include a strong growth of the amplitude of the perturbations. 
Had this occurred the perturbations would have been observed and hence our constraint is 
equivalent to ensuring that such instabilities occur at scales that are not (yet) observable. 

We showed that the presence of corrections to GR does not affect the evolution of 
tensor mode perturbations in a flat radiation dominated universe, while if only corrections 
are present, there is also no effect on the evolution, during the inflation era (although one 
can expect an affect on their production). Requiring that corrections be small for a specific 
mode and a particular scale (within the radiation era) hence constrains the coefficient of 
R'^''Rfj,u corrections and breaks the degeneracy present in the background equations. 

By considering modes that have contributed to the measured (scalar) power spectrum^ 
we can estimate the constraint on R^^'^R^y corrections, finding it to be improved by factor 
of 10'' — 10^ over current best restrictions. Thus, if inflation is to have occurred in the usual 
manner, we require 

a - 3/3 < C (10^°) , a<C(lO^^) , (6.1) 

which is a remarkable improvement on current constraints which give a, /3 <^ O (10™). 

Finally it is important to note that the presence of ghost in theories with such 4**^ order 
corrections is well known 58|] and clearly indicates that the theory cannot be a fundamental 
theory (and certainly cannot be quantised in any standard way). Thus one should consider 
such corrections as an effective theory, which is valid in some range of scales and replaced 
by some more complete theory outside that range. For example, such corrections occur 
in the asymptotic expansion of certain Non-commutative Geometry theories [55|, (in that 
case a = 3/3) and beyond a specific scale the space-time is no longer even approximately 
commutative. Here we have shown that the presence of ghosts is not felt at the level 
of background cosmology, however it is crucial for the evolution of perturbations and the 
constraints produced here essentially ensure that the presence of ghosts would not effect 
observable scales. One expects that higher (i.e. 6**^ and higher 83|, |8J]) order corrections 
to (11.41) would become significant some scale. The presence of ghosts indicate the scale at 
which these higher order terms would need to become significant and essentially mark the 
scale at which the (4*^^ order) effective theory begins to break down. Thus the constraints 
produced here can be viewed as determining the scale at which (11. 4p is a valid effective theory. 
It is important to note however, that here we have assumed a perturbative expansion and 



^ To avoid confusion we emphasise again, the tensor modes have not been measured. We consider these 
modes only because their scalar counterparts have been observed and quantified and if significant en- 
hancement of the tensor modes had occur they would have contributed to the observed power spectrum 
at these scales. 
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it is possible that non-perturbative effects might be significant. 
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Appendix A 

In this appendix we derive the time-time component and trace of the tensor H^^^ given 
in (11. 5p . for the background FRW metric, (12.11) . First note that the non-zero Christoffel 
symbols associated with this metric are, 

f °o = ^ , f}o = ns^, , f °. = n^,, , f f^. = (3)f , (A2) 

where ^^''T^j is the three dimensional Christoffel symbol calculated from the spatial metric 
'jij. In order to evaluate the components of (11.51) we need, 

;p _ ( ffJ. _|_ f a _ f M _ f a \ 

\P ~ ^2 V '^O.O ^ ^ OaJ uO ^ uOJ aO ^ OOJ ua) 

(Jj 

(-pfJ- fa _ -pa _ -pa ffj. \ ( 

where we have defined f^^,^ = R^ya-, which has non-zero components, 

— (OUU' u"\ 
J 00 — {"^^^ ^ ^ ) ^ 

f,o = -^^•(^(^' + 2(^^ + i^)))' . (A4) 



Thus we find, 



Ry^'P^ = 1 [n'" - 2U' {K + Sn^)) . (A5) 
In addition we need ^'°q and R'^^ which are, 

i^'Og = ^ {-en^n' - mn" + {n^ + k) - 2'h'k + n'") , 

R'p^ = ^ [n'" - en^n' - 2v,k) . (A6) 
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Using the expressions given in (12.41) one directly finds the remaining terms in if°Q and 
H^p, eventually giving, 



H 



6 (a - 3/3) 







-2H {H" + 2HH') + m'^ + 2H^K + {H'f - K 



2^^" {n' + K) = P-^ (A7) 



a2 ^ ■ ' 2 

-r^-^{H' + H' + K)=^{p-W) , (A8) 
a"^ 12 

which are the equations given in (I2.5p . 



Appendix B 

Here we calculate the perturbations to -R'^^p, R^-y and R^^.'^'p- In fact, since R'^-p is a 
scalar and we are considering only tensor perturbations, we know that its perturbation will 
vanish ^, however for pedagogical reason we explicitly calculate it here. 

Expanding the covariant derivatives in terms of (derivatives of) the metric we have, 

RP.^p = g"^{d^d0R-T''^^d,R) , (B9) 

where T'^pa is the Christoffel symbol. Perturbing this gives, 

5{R'P.^p) = 5g^^{d^dpR-V\^d^R) 

(d^dpSR - i5Ty^^ d,R - T\^d^5R) . (BIO) 

To calculate this note that 5R = and the non-zero components of the background Christof- 
fel symbols are given in (IA2p and the perturbed versions are. 



(Bll) 



Using this and the expressions in (12. 4p and (13.31) one finds 



D/ 

^ E'm^jR' - {2nEij + R' , (B12) 



— a 



which, as expected, vanishes because our perturbations are traceless i.e. E'^^'jij = 0. 
One can perform a similar expansion of the covariant derivatives in R'^^ to find. 



i?^';, = (/°i?,„)_, + n^/"i?,„. (B13) 



Similarly i?'°.Q, R'^.^, R^o'^p and R^^'^p will vanish as these are either scalars or vectors. 
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Perturbing this one finds (here we only exphcitly calculate the non-zero components), 

~1 — / / TTl/N* 



-R' (E') 



(B14) 

where we have used ( IBIOI) and (13.31) . This is the expression given in (13. 7p . 

The final quantity that we want to perturb is R^y'^'.p- Again, expanding the covariant 
derivatives one finds, 



where 



-f A* — pM pA* _|_ -pt^ pQ po p/i 

J v5 — ^ v:5 — u,5 ^ 5a^ v ^ a 



Perturbing this one finds that the non-zero components of /'^^ and are 



/37 



foj - 



jk — <^ j\k 



where 



= ^ (- iE'% - 2n {EX + {A-K) i?,,) . 
Perturbing (IBlSp we find 

where the (components of the) three tensors Ti, T2 and T3 are given by, 

/ry-i \i — S.„km /-pi rO pO fi pO fi \ 

K-l-l) j — <Jy koJ jm ^ jkJ Om mkJ jO) > 

= in' -n^ -K) e\ , 



(Ts)^ 



mO 



^ jk^J Om ^ jk^J nm "-^ mkJ jO mk^ J jO mk'^ J jn 



(B15) 

(BI6) 



(BIT) 
(B18) 

(B19) 



(B20) 



(B21) 



(B22) 
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Putting these all together one finds 

6 (ir.p^ = In [w -n^-K) {E'y^ + 1 {-o^ - 2ndt + 2n^ + a) , (b23) 

which is the expression given in { \3.7\i . 
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